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Abstract 

We consider well-balanced numerical schemes for the following ID scalar conservation 
law with source term: 

dtu + d^f (n) + z' (x) b (u) = 0. 

More precisely, we are interested in the numerical approximation of the initial boundary 
value problem for this equation. While our main concern is a convergence result, we also 
have to extend Otto's notion of entropy solutions to conservation laws with a source term. 
To obtain uniqueness, we show that a generalization, the so-called entropy process solution 
(see 1 7 1), is unique and coincides with the entropy solution. 

If the initial and boundary data are in L°°, we can establish convergence to the en- 
tropy solution. Showing that the numerical solutions are bounded we can extract a weak*- 
convergent subsequence. Identifying its Umit as an entropy process solution requires some 
effort as we cannot use Kruzkov-type entropy pairs here. We restrict ourselves to the 
Engquist-Osher flux and identify the numerical entropy flux for an arbitrary entropy pair. 
By the uniqueness result, the scheme then approximates the entropy solution and a result 
by Vovelle then guarantees that the convergence is strong in L^, 1 < p < co. 
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In this paper we will consider well-balanced numerical schemes for the initial boundary value 
problem 



(1) 



where Q := ]xi,Xr[, f^r := i7 x ]0, T[ and 

uoeL'^iQ), Ur,uieL'^i]0,T[), feC\R,R), 

are given data. Furthermore we define 
and assume that 

D G (R) , D (R) = R, inf > 0. (4) 

IR 

For the corresponding initial value problem Greenberg et al. [9| have developed and investi- 
gated a well-balanced scheme which is much more efficient than standard schemes. In particu- 
lar stationary solutions can be approximated by these schemes with less numerical investment 
than by the classical schemes. Further results including error estimates and numerical exper- 
iments for the initial value problem have been obtained by Gosse |8|, Perthame [3|. Corre- 
sponding results in multiple space dimensions can be found in f2J and [4J. 

If 6 = there are many results for the initial boundary value problem concerning existence of 
entropy solutions (see Q, O) and convergence of numerical solutions to the entropy solution 
ifTSl . Now in this paper we want to analyse the combination of both, i.e. well-balanced schemes 
for the case 6 7^ 0. 

The outline of this paper is as follows. In Section |2l we will repeat the definition of entropy 
and entropy process solutions. The numerical scheme will be described in Section|3l The main 
result and its proof will be given in SectionlH Finally in Section|5]we will show the advantages 
of these well-balanced schemes in some numerical experiments. 



2 Entropy Process Solutions 

In lfT4l Otto has defined boundary entropy pairs for problems without source terms, which we 
will generalize here to problems with source terms. As in fT?] we will define boundary entropy 
pairs as follows. 
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Definition 2.1 (Boundary Entropy Pairs) 1. Let rj G C^(]R) be convex, q G C^(]R) and 
q' = rj'f. Then [i], q) is called an entropy pair for the partial differential equation in (0). 

2. Let H G C^(]R^), Q G C^{B?), {H{-, w), Q{-, w)) be an entropy pair for all w e R so 
that 

H{w,w) = diH{w,w) = Q{w,w) = Vw G R 

Then (H, Q) is called a boundary entropy pair in the sense of Otto for the partial differ- 
ential equation in (Q). 

3. Let (?7, q) be an entropy in the sense of(l) and w G R such that 

T] (w) = T]' (w) = q (w) = 0. 
Then {rj, q) is called a boundary entropy pair for the partial differential equation in 0. 

There is a close and simple relation between the boundary entropy pairs in the sense of (2) and 
(3) in Definition im 



Remark 2.2 L If{H, Q) is a boundary entropy pair in the sense of Otto (see Definition \2.1\ 2)) 
and w G R then 

7] (u) := H {u,w) , q{u) := Q {u,w) 
form a boundary entropy pair (see Deftnition \2.1\ 3 ) ). 

2. If {rj, q) is a boundary entropy pair (see Deiinition \2.1\ 3)) and w G R such that rjiw) = 
rj'iw) = q{w) = then 

H (u, v) := Tjiu — iy — w)) , Q {u, v) := q{u — {v — w)) 

form a boundary entropy pair in the sense of Otto (see Deiinition \2.1\ 2)). 

Now we use the boundary entropy pairs from Definition |Tlj3) to define an entropy solution of 
the initial boundary value problem ([T]). 

Definition 2.3 (Entropy Solution) Let u G L^iyLx), C := ||u||L°°(nT) and for all boundary 
entropy pairs {rj, q) and all ip G C^iVt x [0, T[), > 0, let 



T 




rj (u) dtip + q (m) dx^p — rj' (u) z' (x) b (n) Lp dx dt 

+ / rj{uo{x))ip{x,0)dx (5) 



where 



Lip (/) / V (ur (t)) ip {xr, t) +rj (ui (t)) (p {xi,t) dt > 0, 
'-c,c] Jo 

J. \f{u)-f{v)\ 

Lip (/) := sup 



[-C,C] -C<u,v<C \U — V\ 

Then u is called an entropy solution o/(E]). Ifz' = this definition corresponds to the definition 
of entropy solution by Otto (see IJ2J}. 
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The basic idea of the convergence proof for the numerical solutions consists in replacing the 
exact solution m in (|5l) by the numerical solution uax as defined in Definition l3.ll First we will 
prove that uax remains bounded for Ax and then we have to control the limit in (jSj) for 
Ax 0. This can be easily done on the basis of the following lemma. 

Lemma 2.4 (Eymard, Gallouet, Herbin ^) Let {uk)ke¥i C. L°°{Qt) be a bounded sequence. 
Then there exists a subsequence {uk')k'e¥S anda function /i G L°°(firx]0, 1[) such that for any 
g e C°(R) we have 

g {uk') / g {fi {x,t, a)) da weak*. 
Jo 

This gives us the motivation for the following definition of the entropy process solution. 

Definition 2.5 (Entropy Process Solution, see [7]) Let /i = fi{x,t,a) be an L°°{Qt>^]0,^[)- 
function, C := ||Ai||L°°(nTx]o,i[) <^nd for all boundary entropy pairs {rj,q) and all ip G C^(f2 x 
[0,T[), > 0, to 

1 

Tj {ji) dtif + q (/i) — rl (fi) z' (x) b (/i) ip da dx dt 



Jn Jo 



+ l^viuoix))(p{x,0)dx (6) 

-c,c] JO 

Then fi is called an entropy process solution o/(|7]). 



+ Lip (/) r r]{ur (t)) ip{xr,t)+ri{ui (t)) if {xi,t) dt>0. 
\-r;r] Jo 



Now we are going to prove that the entropy process solution is unique. Later in Section H] 
we will show that the numerical solutions converge to an entropy process solution. Then the 
uniqueness will give us the existence of an entropy solution. 

Theorem 2.6 (Uniqueness of the Entropy Process Solution) Let fi,u G L°°(r2T'x]0, 1[) be two 
entropy process solutions o/dTJ) with respect to the inital data uq G L°°{fl) and to the boundary 
data ui, Ur G L°°(]0, T[). Then = z/ a.e. on firx]0, 1[ and 

u(x,t):= / fi(x,t,a) da 
Jo 

is an entropy solution o/dU). 



Notice that this theorem also proves uniqueness of the entropy solution. 

In order to prepare the proof of this theorem we need the following definition and some lem- 
mata. 
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Definition 2.7 ( Semi-Kruzkov Entropy Pairs) For A; G R define the semi-Kruzkov entropy pairs 

Tit (.) := {s - kf , g+ (.) := sign+ (s - k) (/ (s) - / (A;)) 
% (s) := - , Qk (s) ■■= sign" (s - A;) (/ (s) - / (A;)) 
where := sign'^(s)s, s" := sign~(s)s and 

1 for s>0, . - . X _ f for s>0, 
/or s < 0, ^^^^ ~ I -1 /or s < 0. 

It can easily be shown that for any semi-Kruzkov entropy pair {rj, q) and any u E [— C, C] we 
have 

\q{u)\ < Lip {f)v{u) . 

[-C,C] 

Now it turns out that we can use the semi-Kruzkov entropy pairs instead of the boundary en- 
tropy pairs in ©. This will be made precise in the next lemma. 

Lemma 2.8 Let fi G L°°(i7Tx]0, 1[), fi = fi{x,t,a) be an entropy process solution o/(|7]), 
C := ||/i||L°°{QTx]o,i[)' k E H and {j]t,qt)> iVk^Qk) semi-Kruzkov entropy pairs in the 
sense of Deiinition \2.7\ Then we can replace {rj, q) in ® by {rj^ , qt) (^nd {rj^ , q^) respectively, 
i.e. for all G C^{VL x [0, T[), (y9 > 0, we have 

III (/^)5t¥' + g^ (/i)5x¥'dadxdt 
Jo JnJo 

— / / sign^ (/i — A;) z' (x) b (fi) (f da dx dt 
Jo Jn Jo f^-. 

f + 

+ / Vk {uo{x))ip{x,0) dx 
Jn 

+ Lip (/) / 7]t (Ur (t)) ip (Xr, t) + 7]^ (Ui (t)) if {xi,t) dt > 0. 

[-c,c] Jo 

Proof Let pk^s (x) := {x — k)'^ + {x — k)^ and define 

r x<k 

vts (^) •= \ Pk,s (x) k < X < k + 6 , 

qtA^) ■■= r%V de 

J k 

Then (r^^^, q^s) ^ boundary entropy pair and we have uniformly in x 

\VkJ i^) - sign"^ {x-k)\ -^0, 

Since {i]ts, Qks) ^ boundary entropy pair we can put it into Q. For 5 — > we obtain (jVl). In 
a similar way we can prove Q for (% , ). 



sign^ (s) 
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The following lemma concerns the trace of an entropy process solution and follows the ideas 
of Otto (see lIH). 

Lemma2.9 Let ji G L°°(firx]0, 1[), jj, = fi{x,t,a), be an entropy process solution o/(|7]). 
Then for all V E L°°(]0,T[) and for all ^ E L^{]Q,T[), (3 >Oa.e., we have: 

ess lim / / sign^ {fx — v) (/ (/i) — f (v)) P da dt 

X^Xr Jo Jo 

>- Lip (/) r {uAt)-vfpdt, 

l-c,c] Jo 

ess lim / / sign^ (/x — v) (/ (//) — / (f )) /5 da dt 

x[xi Jo Jo 

< Lip (/) r {ui{t)-vfpdt. 

[-c,c\ JO 
The essential limits exist. 

Proof Let w G Q be fixed, (77+, g+) be a semi-Kruzkov entropy pair and let /3 G C^(]0, T[), 
/? > 0, and G C^i^), ip>0. Then we obtain with Lemma 



/ (y") P (i) da dt^"' (x) dx 
m Jo Jo 

< r f f (/i) it) ^ (x) 

JO JnJo (8) 

— sign"^ (/i — w) z' (x) h {ji) [3 (t) tp (x) da dx dt 
<C [ ^p{x) dx<-C [ (x) dx 



where C = C{^, z\ 6, (3, w, T). But this implies the existence of a set E' of measure zero such 
that 



<lwil^)P (t) da dt -Cx (x G Q) (9) 

Jo Jo 

is non-increasing in i7 \ E'. Furthermore there exist a set E" of measure zero which only 
depends on n and w, so that the function in Q is bounded on i7 \ E". Therefore the function 
in ® is monotone and bounded on 1] \ {E' U E") and the essential limit exists, i.e. 



lim r [' q+ (/i) /3 (t) da dt = lim C q+ (/i) (3 (t) da dt. (10) 

xr Jo Jo r^f^ „ Jo Jo 



ess 



Similar as in ^ we get for G C^{]xi, x,.]) 

-Iff 

Jn Jo Jo 



x^E'uE" 



(Iw P (t) da dt^' (x) dx 

T (11) 



< C / ^(x) dx+ Lip (/) / r]+{ur{t))(3{t) dt^(x, 
Jn [-c,c] JO 



By convolution with a nonnegative kernel we can also apply dTTT) to i^eix) := ^{x — (xr — e))' 
e > 0. Therefore we obtain 



ess lim (/u) (3 (t) da dt 

X^Xr Jo Jo 



- lim - r r [' g+ (/i) P (t) da dt dx (12) 

eiO e Jx^~-e Jo Jo 



eiO e Jxr~e Jo Jo 

< Lip (/) rrj^{ur{t))P{t) dt 



-C,C] 







By further approximation we can also get (fT2h for /5 G L^(]0,T[), /3 > a.e. Now let v E 
L°°(]0, r[), such that v has only a finite number of values Wi e Q, i.e. 

AT 

v = ^wdA,, AnAj = (!} (t^j), U /l,=]o,i[. (13) 

1=1 i=l,...,N 

Let /5 G L^QO, T[) and /3j := lAiP- Now using and Pi in dT^ and summing over i we obtain: 

ess lim r [' sign± (fi - v (t)) (/ (fi) ~f{v{t)))P (t) d« dt 

x]Xr Jo Jo 

.T (14) 

>- Lip (/) / {Ur{t)-V{t)f P{t)di. 

[-c,c] Jo 

Since any -function v is the limit of such elementary functions of type (fT?l we obtain (fT?b 
also fort; G L°°(]0,T[). 

For the semi-Kruzkov entropy pairs {r]~, q~) and the limits x | we proceed analogously. 

Now we have to show that the entropy process solution has the correct initial data. Again the 
proof is similar to the one by Otto without source terms (see lIT^ . Chapter 2, Lemma 7.41). 

Lemma 2.10 Let /x G L°°(firx]0, 1[), /x = fi{x,t,a), be an entropy process solution o/(|7]). 
Then we have: 

ess lim / / |yU (x, t, a) — (x) | da dx = 0. 



no Jii Jo 

Proof Let ^ G C^(fi), choose (77+, g+) and (r^^, g^) in dTJ) and for e > use 



{x,t) :-- 



1 ■ 
1 - -t 

e . 



l[o,.] {t)tp{x) 



as test function. Adding those two inequalities we obtain 



e JO Jn Jo 
1 r f 



/ / / — \fi — w\ijj (x) + F {n,w) {e — t) ^jj' (x) da dxdt 
Jo Jn Jo 

sign (yU — w) z' (x) b (/i) {e — t)ip (x) da dx dt 
+ / |no (x) — w\iIj (x) dx > 0. 



e Jo Jn Jo 



where 

F(a,6) :=sign(a-6)(/(a)-/(fe)) Va, 6 G R (15) 
Now we take the limes inferior with respect to e and obtain 

esslimsup / / \fi — w\ijj (x) dadx < / \uo (x) — w\ ^p (x) dx. (16) 
ij_o Jfi JO Jn 

Similar as in (fT4l) we get that (O for all w G R implies 

esslimsup / / |/i — f | ^/^ (x) da dx < / \uo (x) — v\ (x) dx (17) 
no J^ Jo Jn 

for all V G L°°{n) and all e L^{n), > a.e. 

Choosing uq for v in (fTTl and using ^ = 1 as test function we get the statement of the lemma. 

The proofs of the following two lemmata are similar to those in ifTSi and therefore we omit 
them. They are mainly based on Lemma l2!9l 



Lemma 2.11 Let ji G L°°(i7Tx]0, 1[) be an entropy process solution o/(|7]). Using F given by 
(1751) we have for all ^ G C^{rtx]Q, T[), ^ > 0: 

/ / / \fi — k\dt(p + F {fi,k) da dxdt 
Jo Jn Jo 

— / / sign (yU — k) z' (x) b (/i) ip da dx dt 
Jo Jn Jo 

+ esslim / / F {fi,Ur (t)) (p {xr,t) dadt 

X^Xr Jo Jo 

— esslim / / F {fi,ui (t)) (p {xi,t) dadt 

xlxi Jq Jq 

- r F {ur (t) , k) (p {xr, t) dt+ r F {ui (t) , k) LP {xi,t) dt > 0, 
Jo Jo 

where /i = t, a) and ip = ip>{x, t). 

Lemma 2.12 Let /i, G L^iVtT^]^, 1[) be two entropy process solutions o/dU with respect 
to the initial data Uq G L°^(fi) and the boundary data ui,Ur G L°°{]0,T[). Then donoting 
fi = fi{x, t, a), V = t, (3) we have for all nonnegative (p> G C^(fix]0, T[), (p> = {p{x, t): 

/ / / \fi — iy\ dt^p + F u) dx^p d(3 dad{x,t) 

Jnr Jo Jo ,io\ 
1 rl ^^^^ 

sign (yu — u) z' (x) {b (/i) — b (z/)) (p d/5 da d(a;, t) > 0, 



inr JO Jo 
where F is given by f l75l) . 
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Now we are ready to prove the uniqueness of the entropy process soiution. The proof without 
the source term can be found in ifTSll . Theorem 2. For controlling the source term we use similar 
ideas as in [1 J, Theorem 5. 

Proof 2.13 (Proof of Theorem\Z^ Let G C^{]0,T[), ^/^ > 0. In (HH) we choose ^ 
tp{t) as a test function and obtain with /i = t, a), v = z/(x, t, (3): 

1 rl 

jji — i'\ dtip if) d[3 da d(x, t) 

(19) 

sign (/i - z/) z' (x) (6 (/i) - h (z/)) (t) d/5 da d(x, t) > 0. 



Q.T JO Jo 



Qt jo Jo 



Using 

g(t):=[ [ [ \n{x,t,a) - iy{x,t,(3)\ dpdadx fortG]0,T[. 
Jn Jo Jo 

([T9b can be written as 



g{t)dti){t) dt< Ik'ILoc(f^) ||&'||ioo(i^) g{t)i){t) dt. 

Integration by parts on the right hand side implies that there exists a set E c]0, T[ of measure 
zero, such that 

t^g{t)- \\z'\\l^^^^ II&'IIl-(r) / 9 (0 
is non-increasing on ]0, T[\E, i.e. for all ti, t2 g]0, ti < t2, we have: 



(t2) < g (ti) + lk'|Loc(f,) ||&'ILoc(i,) / % (0 (20) 
Now we can apply the Gronwall Lemma and (l20b implies 

||/i(-,t2,-) - ^^(-,^2, ■)|lLi(nx]0,l[2) 



< ||;,(.,ti,.)-z/(.,t„.)||^i(^,]o,ip)ell^'"-°°(")"'''"--(-)^*^-*^^ 



(21) 



Since the entropy process solution respects the initial data in the sense of Lemma l2.10l we have: 

ess limsup (-, h, ■) - z/ (-, ti, OH Lif^xio ip) 
hio V J ' L ; 

< essUm -'«o|lLi(f^x]o,i[) + 11^ ^i' O - Mo|lLi(nx]o,i[). 

= 0, 

and therefore from (1211) for almost all t2 G ]0, T[: 

||/i(-,t2,-) -'^(■,^2,-)llLi(nx]o,i[2) <0- 
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But this proves fi (x, t,a) = u (x, t, 13) for almost all (x, t) G VLt and almost all a, {3 G ]0, 1[ . 
Therefore fi and u do not depend on a and /? respectively. Setting 

u{x,t):= / yu(x, t,a)dQ; 

JO 

we get for almost all (x, t, a) G n^xjO, 1[: 

/i (x, t, a) = M (x, t) = u (x, t, a) . 
Then, by Definitions 12 . 31 and 12 . 51 we see that u is an entropy solution of ([T])- 



3 The Weil-Balanced Scheme 



In this section we will describe the well-balanced scheme for the initial boundary value problem 
([T]). It was originally developed by Greenberg et al. in |l9| for the initial value problem. For 
the Engquist-Osher numerical flux Perthame et al. (see |3|) could prove convergence of the 
numerical solution to the entropy solution. Here we are going to generalize this result to the 
initial boundary value problem ([T]). 

First let us fix the notation for the discretization. Let Ax > such that G N, =\xi, Xr[, 
ji G Z and Xj := x/ + (j —ji + 1) Ax, Cj :=]Xj_ i , x^j^i [, choose J C Z such that VL = Ujg j Cj 
and define jV- := max J. 



O// Jl J Jr T 

Let At > such that A^^ := ^ g N and t" := riAt. For any function ip G C^{B?) we define 
V?" := ixj,t") and the piecewise constant function 

^(x,t) ■.= ^] te [t",t"+'), xeCj. 

The main idea for the well-balanced schemes consists in the following fact. If t> is a stationary 
solution of (dJ) then v satifies 

{D {v {x)) + z (x)) = for a.a. x G R, 

(for the definition of D see Q) which is equivalent to 

D{v) + z{x) = c eR for aU x G R. (22) 

This property is the main building block for the numerical scheme. 
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The standard form of a numerical scheme in conservation form for the partial differential equa- 
tion in dl]) is 



At 



g[u],„u]))-Atz'^u] for jG J (23) 
where g is a numerical flux and Zj a discretisation of z' in Cj, e.g. the average of z' on the cell. 



This scheme is very inefficient, especially for the approximation of stationary solutions. The 
well-balanced schemes are much better and are defined as follows. 



Definition 3.1 (Weil-Balanced Scheme) Let 



and define for j ^ J by 



u 



u. 



Similarly let Zj be given by 



1 
1 

r = ^ ^rit) dt for J > jr. 

At Jt- 

Zji for j < ji, 

— z{x) dx for J e J, 
Ax JCj 

Zj, for j > jr. 



(24) 



(25) 



Now assume that {u^)j<z^ is already defined. Then the values m"^ for the new time step of the 
well-balanced scheme are given by 



At 



— (g {ul u]^,^_) - g {u]^,^^, u])) for j e J 



where, due to (EU, u^+i- "^^^ defined by 



D + Zj = D (uj^i) + Zj_i, 

D (mj+i -) + Zj = D {uj+i) + Zj+i. 



(26) 



(27) 



Using the discrete data {u'^)j£w.fi<n<NT define the numerical solution G L°° (R x [0, T[) 
by 

Uax {x,t) := u] for t E [t", T+M. n < Nt, xE Cj, j E Z. (28) 

Remark 3.2 Under condition (|?J there always exist unique solutions and Wj+i,- of X27\) . 
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Remark 3.3 If the values are "stationary", i.e. D{uJ) + zj = D{u'j_i) + zj^ifor all j G J, 
then we obtain 

as solutions of^Tl and therefore u""*"^ = for all j G J. 

Now we assume that the numerical flux is given by the Engquist-Osher flux 

g (m, v) = r /'+ (0 - r (0 d^ + f (0) for u, v G R. (29) 

JO JO 

and set 

H (u, V, w) := V — - — ig (f , w) — g (u, v)) for u,v,w G R. (30) 
Ax 

Before proving the boundedness of the numerical solutions we will prove that the numerical 
values obtained by the Engquist-Osher scheme remain in the convex hull of its data. 

As in [3 1 we use the following "density-function" in order to move the integration bounds to 
the integrand: 

1 for < ^ < s, 

Xs iO :=<!-! for s < e < 0, 
for otherwise. 



Now we can write the Engquist-Osher scheme (l30b as an integral by defining 

At 
Ax 



h (u, V, w) (0 := Xv (0 - ^ (/'^ (0 Xv (0 - r (0 X. (0 



(/'^ (0 xu (0 + r (0 X. 



(31) 



Ax 

for all u, v,w,^ G R. Then (l30b can be written as 

H{u,v,w)= / /i (m, f , w) (^) for M, w G R. 

This form can be used to show that the values produced by the Engquist-Osher scheme are the 
convex hull of its arguments. 

Lemma 3.4 Let h : R? ^ (R) be defined by K = [—C, C], u,v,w ^ K and assume 
that the CFL-condition ^ 

holds. Then for all C, G H we have 

h {u, V, w) (0 G conv {xu (0 > Xv (0 > Xw (0} > 
< sign (^) /i (m, f , w) (^) < 1 

anJ consequently 

|if (m, f , w)| < max , |t>| , |w|} . (32) 
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Proof The proof uses only standard means, see also Q. 



In the next step we will prove the boundedness of the numerical solutions. 

Lemma 3.5 (Boundedness of the Numerical Solutions) Let (|2l), (0), f l29l) be satisfied, 



M := max 



L°°(n) ' ll^dlL°°(]0,T[) ' II'"'-IIl°°{]o,t[)} ' 



j\^g2r||f)'||^cx)(]R^) Ik' II LOO (Q) 



^ ^^ 11^ llL°°(n) ^4r||fe'||,^,„,lk'||,,oo,o, 
.+i)\W\\ 



32T(||6'|| 
-(Ax /~ih.x 



+ \h (0)1 e iii°°(iR.)^-'Jii" iii°°(n). 



anJ we assume the CFL-condition Lipj^A:r (/) < 1- Then we obtain for all n At 



sup \u^A < 



Ax 



Furthermore we have for all {n + 1) At < T: 



sup |<_i , I < sup Im", 1 _| < C^^. 



Proof Let m, t>, w G IR be arbitrary. Then we have with bo := 6(0): 



l^(M,t^) -^(m,w)| 



/ f " (^) ds < / /' 



ds 



</ |f(s)|ds=/ |D'(s)6(.)|ds 
<ll&'lloo/ D'{s)\s\ds+\bo\ I D'{s)ds 

J[v,'w\ J\v,w\ 

< [||6'|Lmax{|t;|,|«;|} + |6o|] / D' [s) 

< [\\h'\\^ms.^{\v\Aw\} + M]\D{v)-D{w)\. 



The mean value theorem implies: 



< 



and therefore: 



< 



'J+1 



Aa;- 



infRD' 
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l&ol 



Aa; II 2;' I 



In a similar way we obtain: 



u 



'j-i 



Ax\\z'\ 



Now by induction with respect to n G N such that nAt < T we will show: 



sup 



OO II 1 1 oo 



+ 



|^^|g2nA<||6'|U+l)||.'|| 



(38) 



which will prove the statement (l35l) of the lemma. Obviously (l38t holds for the initial data, i.e. 
n = 0. 

Assume that (l38t holds for n G N, (n + 1) At < T. Since the maximum principle holds for the 
Engquist-Osher flux g{u, v) (see (l32t ). the CFL-condition holds and Lip^A^: {g) < Lip^^A^ (/) 
we obtain for all j E J: 





< 







Ax\\z' 



u 

+ 
+ 



< sup 



+ 2At 



sup 



Mo 



+ 









^1 




Kl 



< max < sup 



+2At|6o||k'IL + 2AxAt||6'||^ 



I oo II II oo/ 



If max{supj-gj = we have 

< 



u 



n+l 



+ - 



L°°(]0,T[) 

Ax lU' 



g2At||6'|U||y|Ug2nAt||fe'||^ll^'„o, 



>1 



infRD' 



<^g4At||6'|UP'||^g4nAt||b'|UP'„oo 



>1 



< Me^("+^)^*"'''l'-ll^'ll- + ll^'IL ^4(n+l)At||fc'| 



CXD I! I! oo 



infR^' 



+ 



|5g|e2("+i)^*(ll^'IL+i)P'IL. 
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In a similar way the case maxjsupjgj \uf\, = can be handled. 
If maxjsupjgj \uf\, = sup^gj we use (l38t to obtain: 



< Me2"^*ll'''ll-ll"'ll-(l + 2At||6'| 



Ib'll 1 

CXD II II 00/ 



^ ^^\W\\oo ^AnAt\\b'\U\z'\\^ 2At ||/|l ) 

inf^I^' °° °° 

+ |6o|e^"Mll^'ll^+i)ll^'IU (l + 2At\\b'\\^\\z'\\J 

/||2 



+2At|6o| ||/||^ + 2AxAt||6'| 



L2 



+ 



|5g|e2{n+l)Ai(||6'|U+l)ll^'lloo_ 



Therefore we have proved (l38t for nAt < T. The estimates in (l36t now follow from dJTt . (l38b 

and (n + 1) At < T. 

Notice that the second term in (l34l) is necessary, although it is missing in the proof for the initial 
value problem in |3|, Lemma 3.1. This is shown in the diplom thesis [ 13|. 

A main tool for the convergence proof in SectionlUwill be a cell entropy inequality. This means 
for an arbitrary entropy pair (77, q) we have to find a numerical entropy flux G, such that for 
C > and for sufficiently small 

T] (H (m, V, w)) -viv) + ^(G (v, w)-G (u, v)) < (39) 

Ax 

holds for all u,v,w E [—C, C]. 

To prove the inequality in (l39l) we will need a lemma by Brenier (see f5\). 

Lemma 3.6 (Brenier /|5j/j Let f E (R) such that < sign (0/(0 — ^fo^ all C, E H and 
let h E G"^ (R) be convex and Lipschitz-continuous. Then we have: 



h (/^ / (0 de) - (0) < h' (0 / (0 de 



Now we can prove the cell entropy inequality (l39t . which generalizes a result due to Perthame 
et al., see [3|. 

Lemma 3.7 Let \u\ , \v\ , \w\ < G E H, K := [— C, G], (77, q) be an entropy pair and assume 
the CFL-condition 

II/'IIl-w|^<1- (40) 
15 



Then G : R x R ^ R given by 

G (w, V) := r V (0 /'+ (0 - r V (0 (0 de + g (0) (41) 
JO Jo 

is a consistent numerical entropy flux, such that the cell entropy inequality d^PI ) is satisfled. 

Proof It is easy to check that G is a consistent numerical entropy flux. Since |u| , |f | , |w| < C 
and due to (l40b \H {u,v,w)\ < C, we can assume without loss of generality that rf is Lipschitz- 
continuous. 

Due to the CFL-condition (l40b and Lemma l34l we have for all ^ G R that 

Q<sign{i)h{u,v,w){i)<l 

and by Lemma IT6l we have 

r/ h (m, V, w) (0 de) - 7] (0) < 1^ r/' (0 (m, i;, w) (0 d^ 

Then using dSTT) and (|4T]l on the right hand side, the cell entropy inequality (l39b follows. 



4 Main Result and Convergence of the Numerical Solutions 

In this section we will present the main result of this paper, the convergence of the numerical 
solutions of the well-balanced scheme (see Definition l3.1l) to the entropy solution, and its proof. 

Theorem 4.1 Let 0, (0) and d^PI) be satisfled, let (max)ax>o be the sequence of numerical 
solutions o/0, in the sense of Definition \3 .l\ and assume the CFL-condition Lip^Ai. (/) < 
1, where C^^ is deflned as in (1341) . Then there exists au & L°° (Qt), such that 

uax u in (^t) for all 1 < p < oo (42) 

and u is the unique entropy solution o/0. 

Proof see Proof l4.5l 

For the proof we need the following lemmata. 

Lemma 4.2 Let (r/, q) be a boundary entropy pair, let w E H such that rj {w) = rj' {w) = 
q {w) = and let G be deflned by X4l\) . Then we have for all u,v E H 

G (n, v) = r (0 /'+ (0 de - r/' (0 (0 de (43) 

Jw J w 

If additionally u.,v,w G [— C, C], C > 0, holds, the following inequality holds true: 

- Lip {f)v{v)<G{u,v)< Lip {f)r^{u). (44) 

[-C,C] [-C,C\ 
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Proof Since q {w) = we have 

g(0)= fq'iO de= f\'{OnO 

= f ri (0 r (0 - f V' (0 /'^ (0 
Taking (|4T1) into account, we obtain (l43b . 
Since t] E (R) is convex and r^' (w) = 0, we have 

sign {u — w)ri' {^) > for all ^ G [u, w] . 
Therefore, as /'^ > 0, we see that 

/" r]' iO iO de = / sign {u - w) (0 /'+ (0 d^ > 0. 

Jw J[w,u] 

On the other hand i] (w) = and we can estimate: 

/\'(0/'^(0 de < ll/1Loc([_c,c]) / sign(n-^)V(0 d^ 

J«) " J [w,u\ 

= Lip 

[-C,C] 

Applying the same arguments we see that 

0< d^< Lip if)r]iv) 

Jw [-C,C] 

and using (l43t we obtain (l44l) . 

Lemma 4.3 Let (0), 0, (1291 ) anJ ce/Z entropy inequality (1391 ) satisfied. For At, Ax > 
assume that < A anJ /e? iu'j)o<n<NT,jeJ ^^e numerical solution o/dTJ) m j'eni'e o/ 
Definition I3.il Le? (r/, g) be a boundary entropy pair, w G R such that rj {w) = rj' [w) = 

q{w) = 0, |w| < C for all j E J, n < Nt and \u^_^_i_\,\u^_i^^\ < C for all j G J, 
n < Nt. Then for all <p G C^{R'^) n C^{U x [0, T[), > 0, we have ' 

[ [ V (uax) dtif + q (uax) d^(p - r]' {uax) b (uax) z' (x) ip dx dt 
Jo Jn 

+ r]{uo{x))ip{x,0) dx+ Lip (/) / r] {ur (t)) ip {xj^,t) dt (45) 
Jn [-c,c] JO ^ 

+ Lip (/) rr]{ui{t))ip{x,,,t) dt>0{Ax). 
l-c,c] Jo 
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Proof Multiplying the cell entropy inequality ( l39b by Ax and summing over n < Nt and 
j G J we obtain 

Nt-1 

n=0 jeJ 
Nt-1 

Nt-1 
n=0 j'eJ 

n=0 j£j 



(46) 



We will now treat the four sums in (I46t seperately. The first term is handled in a standard way 
by an index shifting argument: 



Nt-1 

E {v{ur')-ri{u-))^ 

n=0 

Nt-1 

Nt\ Nt-1 ( OA -1 

n=0 

Nt-1 



T. n (< - + .) («f ) - (.«) (47) 



n=0 



where < At ll^tVlL + 2 \\dM\J ■ 

Using Jensen's inequality we deduce that 

1 f 

and from (l47t we obtain 

7Vt-1 



n=0 



(48) 



n=0 



Now we turn to the second term of Using an index shifting argument and (HH) we 
obtain 

T.{g{^1^1^i)-G{u]_„u]))^] 
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> -Ax 5: g («^") (S.y.); - r, {ul^,) Lip (/) 



-V Lip (/) + S,,^. 

[-C,C] 



where \£x,n\ < Cx,n^x with 



IIVIIlo.([_c,c]) ll/'llL^(hc,c]) (1^1 I^'V'IL + ll^-'^llc 



Again using Jensen's inequality to get 



1 f* 



1 /■* 



we arrive at the following inequality: 



>-AxJ2q {u]) {dxV)] - Lip (/) — / (t)) dt 

- Lip (/) — / 7^{ui{t))y,ldt + Sx,^. 

\-C,C] Jt" 



(49) 



Now we turn our attention to the discretisation of the source term. By the mean value theorem 
there is a C]Vi e [m^+i, -] such that 



= - / (0 (0 de 

j+i 

Also by the mean value theorem and (l?7t there is a E [u'j^i, 

^,+1 -z, = D {u]^,^_) - D [u]^,) = D' 
and since inf]R > we have 



so that 



7/" 



Ax 



(50) 
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Since zj^ — zj^^i = Zj^+i — zj^ = 0, multiplying (l50b by (f"' and summing over j G J we obtain 



where 



i6J 



—Ax y^^rj' \u 



D'{^]) Ax 
'^^ D'iu'i) Ax 



n I c 



(51) 



By the assumptions © 4? is differentiable and we can estimate: 



/ (.n\ f (C;) _ , fn\ f (»i) 



I < C*Ax^ 

J IT 



where C* = C* (C, rj, /, D). Using this result we see that \£+^n\ < c±_„Aa; with 



^" l°°([-C,C]) II/1Il°°([-C,C]) II ,|| 

inij^D' IF lloo 117- noo 



In the same way we treat the other part of the source term and obtain 



(52) 



where \£-^n\ < c±^„Ax. 
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Using (l48t . (l49l) . dSTT) and (|52l) in (l46t we arrive at the following inequality: 

Nt-1 

At Ax Y: E V (<) (dt^); + AtAx E E g (tz,") (S.y.) 

n=0 jeJ n=0 jr'eJ 



J 



0<n<JVy 



D'{u]) Ax D'{u]) Ax 



^1 



(53) 



+ Lip (/) E / V {Ur (t)) + V (ui it)) dt > O (Ax) 



[-C,C] „=o 

It is not hard to see that ( 15^ implies (B31) . 

For the final convergence proof we need the following relation between weak* and strong 
convergence: 

Lemma 4.4 Let Q be a bounded, measurable subset o/R", let {vn)^^^ i>e a bounded sequence 
in L°° (Q) and v G L°° (f2) such that for all continuous functions g E C (R) we have 

9 {vn) g (f ) weak * in L°° (fl) . 

Then the convergence Vn —>■ v is strong in L'^ {Q) for 1 < p < oo. 

Proof see ifTSl . Lemma 10. 



Proof 4.5 (Proof of TheoremKB Let 5 > be fixed. Then for Ax < 5 we have ^ < C^, 
where C^^ is defined by (l34l) . Hence the CFL-condition and Lemma l331 provide for Ax < 5 

\\'^Ax\\L°^(nx]0,T[) — ^T- 

This lemma also provides the boundedness of the discrete equilibrium states, as defined in (l27t . 

Let (?7, q) be a boundary entropy pair and let G R such that rj [w) = i]' (w) = q (w) = 0. By 
Lemma l3Jl there is a consistent numerical entropy flux G, so that the cell entropy inequality 
(l39b is satisfied. 

Let V? e C^(n X [0, T[), v9 > 0. Then ip can be extended to a function in C^{U x [0, T[) n 
Cq(R^), which we will also denote as (p. 

Setting C := max{C^, } all the premises of Lemma l4.3l are fulfilled and we can deduce that 
for < Ax < (5 and for all non-negative test functions (f G C^{U x [0, T[) fl Cq (R^), we have 
151). 



21 



Now Lemma l24l allow s us to find a sequence {Axk)ke¥i c]0, 6] and a G L°°(fiTx]0, 1[) such 
that for all 5f G C(R) we have 

g{uAxk)-^f 5' (/"(■,■,«)) da weak* in (fir) . (54) 

JO 

Since r], q, rj'h G C(R) and dtf^, dx^, z'^p G ^^(^7^) we can pass to the limit in (l45t and obtain 
Thus /i is an entropy process solution of ([T]) and by Theorem l2.6l get a m G L°°{^It) such 
that a) = M a.e. for almost all a g]0, 1[. 

Since the entropy process solution is unique by Theorem l2.61 we conclude that any subsequence 
of {uax)ax>o converges to weak* to u as Ax 0. By Lemma we see that {uax)ax>o 
converges strongly to u in Lp{Q,t). 



5 Numerical Experiments 



In this section we will numerically test the well-balanced scheme (see Definition l3.ll) and com- 
pare it to the standard discretization (l23t . Therefore we take a look at the Burgers-Hopf equa- 
tion 

5tM + 5xy + = (55) 

on :=]0,4[, which means D{s) = s. So © is obviously satisfied. For the source term we 
choose 

cos (vrx) for x G 



3 5 \ 

z(x) = ■l^^^^"'^' — -V2'2j'. (56) 
[ otherwise 

Notice that z' is discontinuous and Uc{x) := c — z(x) is an equilibrium of (l55t . 

Test Case 5.1 Assume the equilibrium initial data M2 := 2 — z{x) and the constant boundary 
data ui = Ur = 2. The exact solution to this problem is given by u{x, t) = U2{x). 

Now we use 40 nodes in space and compare the numerical solutions of both schemes (as defined 
in (13.11) and (l23l) respectively) at t = 3. The result is shown in Figure[l] 

The Tables[l]and|2lshow the L^-error for the numerical solution obtained for different grid sizes 
Ax. They also display the L^-error of the numerical solution with respect to the projection of 
the exact solution to the space of piecewise constant functions (Numerical Error). 

In these computations At is computed such that the CFL-condition is satisfied. The necessary 
estimate on is given by Lemma U31 For the standard discretisation (E^ we have 



< sup + At ll-z'll^^ooi'f^^ sup 



where C^^ is defined as in ( B^ . So our estimate for the well-balanced scheme seems far from 
optimal. 
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2.5 



1.5 



0.5 




.0...<j...^...^...^...^...<^..0 ..0..,<j...<j...^...^..^...^ 



♦ Well-Balanced Scheme 
o Standard Scheme 
exact solution 



0.5 



1 



1.5 



2.5 



3.5 



Figure 1: Well-Balanced and Standard Scheme for Test Case lS.ll at t = 3 



Ax 


At 


-Error 


Num. Error 


CPU-Time 


10-1 


5.7-10-'' 


5.02 - 10-2 





0.62 s 


10-2 


3.5 - 10-*^ 


5.00- 10-3 





8.83 s 


10-3 


7.3 - 10-^ 


5.00 - 10-4 





423.2 s 


10-4 


8 - 10-8 


5.00-10-^ 





42 688 s 



Table 1: Well-Balanced Scheme for Equilibrium Initial Data uq = z (x) 



Ax 



At 



L^-Error Num. Error CPU-Time 



2 


10" 


-2 


4.5 


10-^ 


5.07 


10- 


-2 


4.90 


10- 


-2 


2.94 s 


2 


10- 


~3 


1.3 


10-6 


5.09 


10- 


-3 


4.91 


10" 


~3 


97.6 s 


2 


10- 


-4 


1.6 


10-^ 


5.09 


10- 


-4 


4.92 


10- 


-4 


10 885 s 


2 


10- 


-5 


2 


10-8 


5.06 


10" 


-5 


4.97 


10- 


-5 





Table 2: Standard Scheme for Equilibrium Initial Data 2 — z (x) 
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Weil-Balanced Scheme 
Standard Scheme 
- Standard Scheme with 40000 Nodes in Space 



t = 0.25 



Weil-Balanced Scheme 
Standard Scheme 
- Standard Scheme with 40000 Nodes in Space 




t = 0.75 




Weil-Balanced Scheme 
Standard Scheme 
- Standard Scheme with 40000 Nodes in Space 











• Well-Balanced Scheme 




Standard Scheme 




Standard Scheme with 40000 Nodes in Space 





t = 1.75 t = 2.75 

Figure 2: Well-balanced scheme vs. standard discretization for Test Case l5.2l 



Test Case 5.2 Assume the constant initial data uq = 1 and the constant boundary data ui = 2, 
Ur = 1. Again we use 40 nodes in space to compare the well-balanced scheme as defined in 
(13.11) to the standard discretization (l23t . A reference solution was obtained by the standard 
scheme with 40 000 nodes in space. The result is displayed in Figure |5]2| 



Now we turn to the weakness of the well-balanced scheme. For the Burgers-Hopf equation 
(l55t u{x,t) = is an equilibrium for any choice of z. Thus, if we ask the boundary data 
ui = Ur = there are two equilibria satisfying this condition: —z{x) and 0. 

Test Case 5.3 We assume the constant initial data uq = and the constant boundary data 
ui = Ur = 0. the exact solution is given by u{x,t) = 0. Figure |3l displays the numerical 
solutions obtained by the well-balanced scheme for 40, 400, 4 000 and 40 000 nodes in space at 
t = 2.5. 

Table |3] shows the L^-error for this problem, which is not acceptable even for 4 000 nodes in 
space. 
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♦ Weil-Balanced Scheme with 40 Cells 
Well-Balanced Scheme with 400 Cells 

Well-Balanced Scheme with 4000 Cells 

Well-Balanced Scheme with 40000 Cells 



0.5 1 1.5 2 2.5 3 3.5 



Figure 3: Well-balanced scheme for Test Case 15 .31 at t = 2.5 



Ax 


At 




-Error 


CPU-Time 


10-1 


6.14-10- 


-6 


4.388 - 10-1 


0.90 s 


10-2 


6.14-10" 


-6 


3.164- 10-1 


5.23 s 


10-3 


6.14- 10" 


~6 


2.678 - 10-2 


47.9 s 


10-^ 


6.14-10" 


-6 


8.421 - 10-^ 


740.8 s 



Table 3: Well-Balanced Scheme for Initial Data uq = 
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3.5 



3 - 




2.5 - 



2 



1.5 - 



1 - 



♦ Weil-Balanced Scheme with 40 cells 
o standard Scheme with 40 cells 

Standard Scheme with 4,000 cells 
Standard Scheme with 40,000 cells 



0.5 1 1.5 2 2.5 3 3.5 



Figure 4: Well-Balanced and Standard Scheme for Test Case 15 .41 at t = 3 



Notice that the convergence proof for the well-balanced scheme holds for m = 0, even though 
we divided the differential equation in ([T]) by m in the motivation. 

We also tested both the well-balanced scheme and the standard scheme (l23b for a discontinuous 

z: 

Test Case 5.4 Let z be given by the discontiuous function 



and choose the initial data U2 = 2 — z{x) and the constant boundary data ui = Ur = 2. The 
numerical results in this case are displayed in FigurelH 

Note that while the well-balanced scheme converges to the equilibrium quite fast, the standard 
scheme seems not to converge to the equilibrium at all. 
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otherwise 
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